Non-white frequency noise in spin torque oscillators and its effect on spectral 
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We measure the power spectral density of frequency fluctuations in nanocontact spin torque 
oscillators over time scales up to 50 ms. We use a mixer to convert oscillator signals ranging from 
10 GHz to 40 GHz into a band near 70 MHz before digitizing the time domain waveform. We analyze 
the waveform using both zero crossing time stamps and a sliding Fourier transform, discuss the 
different limitations and advantages of these two methods, and combine them to obtain a frequency 
noise spectrum spanning more than five decades of Fourier frequency /. For devices having a free 
layer consisting of either a single Ni8oFe2o layer or a Co/Ni multilayer we find a frequency noise 
spectrum that is white at large / and varies as 1// at small /. The crossover frequency ranges from 
^ 10^ Hz to ~ 10^ Hz and the 1// component is stronger in the multilayer devices. Through actual 
and simulated spectrum analyzer measurements, we show that 1// frequency noise causes both 
broadening and a change in shape of the oscillator's spectral line as measurement time increases. 
Our results indicate that the long term stability of spin torque oscillators cannot be accurately 
predicted from models based on thermal (white) noise sources. 



INTRODUCTION 

In a spin torque oscillator (STO), a direct current pass- 
ing through a reference magnetic layer becomes spin- 
polarized and transfers angular momentum to a second 
magnetic layer that is excited into steady-state oscilla- 
tion. The oscillating magnetization causes an oscillating 
device resistance, through either the giant magnetore- 
sistance effect or the tunneling magnetoresistance effect, 
which in combination with the bias current generates an 
oscillating voltage as the output signal. Interest in poten- 
tial applications of STOs in integrated microwave circuits 
is driven by their rapid frequency tunability, small size 
(< 100 nm), and compatibility with standard semicon- 
ductor processing techniques. Recent reviews cover both 
the physics [1 and possible applications [2| [3j of STOs 
and other devices based on spin torque effects. 

For any oscillator, noise is both an important figure 
of merit for applications and a useful probe of internal 
physical processes. Previous models of STO noise [4]-[9] 
have considered noise driven by a thermal source having 
a white power spectral density (PSD). Perhaps because 
most experiments on STOs have used a spectrum ana- 
lyzer (SA) to measure signals in the frequency domain, 
most previous theoretical work has focused on how fre- 
quency noise determines the width of the spectral line. 
For purely white frequency noise, the relation is straight- 
forward: a constant PSD 5'wh gives a Lorentzian spec- 
tral line whose full width at half maximum is simply 
Ai^wh = 7r5'wh Uni HI]- The situation is more compli- 
cated when the frequency noise is not white. Colored 
noise in STOs can occur both at high frequencies, due 
to an intrinsic relaxation rate that suppresses rapid fre- 
quency fluctuations [7, 9 , and at low frequencies as we 
demonstrate here. The stability of an oscillator cannot 
be described by a single number such as linewidth when 
its frequency noise is colored, and a measurement of the 



noise spectrum is required to accurately predict oscillator 
performance in specific applications and to test models 
of the physical origin of the noise. 

In the next section we present generic equations for an 
oscillator in the time domain and introduce the various 
PSDs we use here. Then we describe our STO devices 
and our measurement techniques, including the use of a 
mixer to facilitate the measurement of signals well above 
10 GHz. Next we describe two methods (both employing 
standard digital signal processing techniques) for com- 
puting the PSD of frequency fluctuations from the volt- 
age waveform of the STO. While each method has dif- 
ferent bandwidth limitations, when combined they yield 
noise spectra extending over more than five decades of 
Fourier frequency /. We present these combined spec- 
tra for two types of STOs, both of which show 1// fre- 
quency noise below / ^ 1 MHz. Finally, we connect the 
frequency noise with SA measurements. The measured 
linewidth is larger than the value implied by the white 
part of the frequency noise spectrum and it increases 
with measurement time, effects that have been seen in 
semiconductor lasers having 1// frequency noise. The 
line shape also changes, becoming more Gaussian at long 
measurement times, and we discuss how this affects the 
interpretation of STO line shape measurements. 

TIME DOMAIN OSCILLATOR MODEL 

The voltage output of a generic oscillator can be writ- 
ten as 



V{t) = [Vo + e{t)] sin pTr^^o^ + ^t)] , 



(1) 



where e{t) is the deviation from the nominal amplitude 
Vo^ fQ is the nominal frequency, and (j)(t) is the deviation 
from the nominal phase 27ruot. From the total phase 



0{t) = 27Tiyot + 0(t) 



(2) 



we define an instantaneous frequency 

1 de 1 d4> 



(3) 



From this equation, it is clear that phase and frequency 
are equivalent, not independent, ways of representing os- 
cillator fluctuations. 

Oscillator noise is commonly expressed in terms of 
the PSD [27 of V{t), (l)(t), or v{t), which we denote as 
Sv{y) [VVHz], S^{f) [radVHz], and S,{f) [HzVHz]. 
The units of each PSD are given in square brackets and 
/ is Fourier frequency. Note that 6'y(i^), the quantity 
measured by an SA with swept frequency z^, includes 
amplitude noise that does not appear in the other two 
PSDs. These PSDs can be measured in various ways; 
we will describe the methods we use below. We will also 
make use of the Fourier identity 



SM) = fsM), 



(4) 



which follows from the fact that frequency is the time 
derivative of phase. 
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FIG. 1: (Color online) Spin valve layers used for the NiFe 
and Co/Ni STO nanocontact devices. In each case, the layers 
extend laterally over a mesa of several /xm on a side. The 
current that generates the spin torque effect flows through a 
metallic nanocontact to the top layer and returns through the 
bottom Cu layer. 



EXPERIMENTAL METHODS 



Our nanocontact STOs consist of a laterally extended 
spin valve structure and a metallic contact of nominal 
diameter 60 nm to 70 nm. Figure [l] shows the two types 
of structures we used. Both structures have a thick layer 
of CogoFeio that serves as the reference layer. We use 
"NiFe" to label devices whose free layer consists of 5 nm 
of Ni8oFe2o and "Co/Ni" to label devices whose free layer 
consists of a multilayer of Co and Ni. With no applied 
magnetic field, the free layer magnetization of the NiFe 
devices lies in the film plane, whereas that of the Co/Ni 
devices lies perpendicular to the film plane due to inter- 
facial anisotropics intrinsic to the multilayer [12 . The 
NiFe devices are from the same wafer as those used in 
[13] and the Co/Ni devices are from the same wafer as 
those used in fT4 . 

We used microwave probes to contact the devices, 
which were at room temperature. The high-frequency 
STO output, after separation from the bias current by a 
bias tee, passed through an amplifier with a power gain 
of 36 dB before entering an SA. We measured more than 
a dozen devices from four different wafers, at a variety 
of applied magnetic fields and bias currents, and we ob- 
served no n- white frequency noise in all cases. Here we 
present representative data from one NiFe device and 
one Co/Ni device. The NiFe device was measured with 
a magnetic field fioHo = LOT applied at an angle of 80° 
from the film plane and a bias current /b = 12.1 mA, for 
which uq = 13.2 GHz. The Co/Ni device was measured 
with jhqHo = 1.2 T applied at 85° from the film plane 
and /b = 7.5 mA, for which uq = 36.5 GHz. We chose 



these conditions to illustrate how devices with nearly the 
same SA linewidth (?^ 11 MHz in this case) can have sig- 
nificantly different frequency noise spectra. 



Our technique for time domain measurements was mo- 
tivated by a desire to measure STO signals well above 
10 GHz using readily available commercial instruments. 
We used the intermediate frequency (IF) output of the 
mixer in an SA, with a fixed local oscillator (LO) fre- 
quency, to translate the input signal V{t) to an IF sig- 
nal ViF{t) centered at 70 MHz. Ignoring the negligible 
phase noise of the LO, a perfect sine wave V{t) would 
appear at the IF output as a perfect 70 MHz sine wave, 
whereas the frequency or phase fluctuations in an actual 
V{t) appear as corresponding fluctuations in Vw{t). The 
advantage of the IF measurement technique is that STO 
signals for any uq within the range of the SA are trans- 
lated to a common, lower frequency at which digitization 
is straightforward. A disadvantage is that the limited 
IF bandwidth prevents the measurement of signals with 
large linewidths. For the data presented here, the avail- 
able IF band was approximately (70 ± 20) MHz and we 
were limited to STO signals with linewidths < 30 MHz. 
We used an oscilloscope to digitize Vwit) at 10^ samples 
per second and to apply a digital 150 MHz lowpass filter 
to reduce preamplifier and oscilloscope noise. Figure [2] 
shows a portion of the resulting waveform. The available 
oscilloscope memory limited the duration of each filtered 
waveform to < 50 ms. 
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FIG. 2: A short section of the IF waveform ViF(t) for the 
NiFe device. The ampUtude variations seen in these data are 
dominated by amphfier noise. 



DATA ANALYSIS METHODS 

In this section we describe two different methods for 
obtaining Sjy{f) from Vwit). We discuss Hmitations and 
averaging considerations in some detail and we show that 
the two methods agree over their common range of /. 
Beyond the points we highlight here, many textbooks 
and other sources (e.g., [15]) contain details of the digital 
signal processing techniques involved. In the next section 
we present composite spectra obtained by combining the 
two methods in order to cover a broader range of / than 
is possible with either method alone while preserving the 
benefits of averaging. 

The first analysis method is the "zero crossing" 
method. As described in ^, we obtain a value of the os- 
cillator phase each time I^if(^) crosses zero, which yields 
(j){t) as shown in Figs.lsFa) andlsFb) for the NiFe device. 
Because the values of (j) are not equally spaced in time, 
one should in principle estimate the PSD using an algo- 
rithm such as the Lomb periodogram [15] that is suitable 
for a time series with irregular spacing. However, in prac- 
tice the variations in spacing are sufficiently small that we 
find no significant difference between the Lomb PSD and 
that obtained by assuming uniform spacing and applying 
conventional PSD algorithms. Thus for the analysis pre- 
sented here we have replaced the actual time stamps in 
each <j){t) trace with values separated by the mean spac- 
ing for that trace. To reduce scatter in the spectra [15], 
we averaged PSDs computed from half-overlapping seg- 
ments of ^(t), each multiplied by a Hann window, to ob- 
tain S(f){f) shown in Fig.lsFc). We chose segment lengths 
of 1 ms and 10 ms to balance the tradeoff between aver- 
aging and frequency resolution. We found it necessary to 
omit the lowest three frequency bins from S(f){f) to obtain 
results that are independent of segment length; thus for 
10 ms segments 5*0 (/) begins at / = 400 Hz rather than 
(10ms)~^ = 100 Hz. (Spurious effects arise when there is 
a large difference between the initial and final values of (j) 
for a segment. Consider a segment where (j){t) varies lin- 



early: multiplying a line by a window that falls to zero at 
each end will create artificially large Fourier components 
near the inverse segment length.) Finally, we use Eq. Q 
to obtain Sy{f) for the zero crossing method, with the 
result shown in Figure IsF 
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FIG. 3: (Color online) Zero crossing method for computing 
Su{f) from ViF{t). (a) Short section of (j){t) from zero cross- 
ings of Vw{t) for the NiFe device, (b) (/)(t) for the entire 50 ms 
IF waveform, (c) Average S(f){f) computed from segments of 
(/)(t) to reduce scatter in the PSD. (d) Su{f) computed using 
Eq. Q. For both (c) and (d), dotted curves are the aver- 
age PSD from 9 half-overlapping, 10 ms segments (points for 
/ > 5 X 10^ Hz are not shown) and solid curves are the average 
PSD from 99 half-overlapping, 1ms segments. S(f)(f) flattens 
near 10 ~^ rad^/Hz due to the noise floor for the zero crossing 
method set by amplifier noise in our setup [16 , which causes 
the upturn in Su{f) at large /. Dashed lines are visual guides 
indicating 1//^ and 1// power law spectra. 

The rolloff in S^{f) beginning at / ?^ 20 MHz qualita- 
tively resembles the expected rolloff due to the intrinsic 
relaxation rate at which an STO returns to its stable 
precessional orbit after a fluctuation [9]. However, in our 



case this rolloff is due to the hmited bandwidth of the IF 
output (it occurs near 20 MHz in ah our data) and does 
not reflect an intrinsic time scale of the STO. For the 
rest of this paper we show Sy{f) from the zero crossing 
method only for / < 20 MHz. 

The second analysis method is the "sliding DFT" 
method. We compute the discrete Fourier transform 
(DFT) of half-overlapping segments of Vwit) and flt a 
Lorentzian peak to the DFT to determine the center fre- 
quency for each segment [28 . This yields a trace of z/(t), 
as shown in Figs.pFa) and[4Fb) for the NiFe device. The 
shortest segment that gave reliable results was 300 ns, 
which for half-overlapping segments gives a value of v ev- 
ery 150ns. We then compute Sy{f) for the DFT method 
from i/(t), again averaging the PSDs of half-overlapping, 
Hann-windowed segments. Since the relative excursions 
in v(t) are much smaller than those in (j)(t)^ there are 
no spurious effects for long segments and Siy{f) is inde- 
pendent of segment length over the entire range of /. 
We could easily obtain 5*0 (/) for the DFT method using 
Eq. Q, but the features of interest here are more easily 
seen in Sy{f). 

In Fig. [5] we compare the average Sy{f) for 10ms seg- 
ments obtained from the two analysis methods. The two 
methods give nearly identical results over their common 
frequency range. Such agreement between two differ- 
ent routes to the same quantity is evidence that neither 
method is distorted by numerical artifacts, and thus that 
both methods reveal the actual frequency fluctuations of 
the oscillator [29 . Comparing frequency ranges in Fig. pj 
the zero crossing method extends up to / = 20 MHz, 
while the DFT method is limited to / < 3 MHz by the 
minimum DFT segment length mentioned above. On 
the other hand, since the DFT method avoids the spuri- 
ous effects related to segment length, it extends to lower 
frequencies than the zero crossing method (for a given 
segment length). 
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FIG. 4: (Color online) Sliding DFT method for computing 
5'^(/) from Vw{t). (a) Short section of v{t) from DFTs of 
half-overlapping segments of ViF{t) for the NiFe device, (b) 
zy(t) for the entire 50ms IF waveform, (c) Average Su{f) 
computed from segments of vit) to reduce scatter in the PSD. 
Dotted curve is for 9 half-overlapping, 10 ms segments (points 
for / > 5 X 10^ Hz are not shown). Solid curve is for 99 half- 
overlapping, 1 ms segments. Dashed line is a visual guide 
indicating a 1// power law spectrum. 



RESULTS AND DISCUSSION 

As mentioned above, we combine spectra from the zero 
crossing and DFT methods to obtain Sy{f) over a larger 
frequency range than is possible with either method 
alone. As shown in Fig. |6J the zero crossing method 
using 1 ms segments together with the DFT method us- 
ing 10ms segments yields Sj,{f) spanning more than five 
decades in /. Both the NiFe and Co/Ni devices show 
the same overall behavior: Sy{f) is constant at large / 
(excluding the rolloff and noise fioor features mentioned 
above) and varies as approximately 1// at small /. We 
found this same qualitative behavior in all STOs we mea- 
sured. For a given device, we have not found a clear 
dependence of either the white or 1// noise on bias cur- 
rent or applied field, but our measurements to date have 
consisted of a broad survey rather than a search for sys- 
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FIG. 5: (Color online) Comparison of Su{f) from two anal- 
ysis methods for the NiFe device. The curves are those for 
10 ms segments shown in Figures [3|d) andlljc), here shown 
for all / < 20 MHz. The agreement indicates that numerical 
artifacts are likely negligible in both methods. 
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FIG. 6: (Color online) Combined frequency noise from both 
analysis methods for (a) NiFe device and (b) Co/Ni device. 
Segment length is 1 ms for the zero crossing method and 10 ms 
for the sliding DFT method. Dashed lines are visual guides 
indicating a 1// power law spectrum. 



tematic trends. One clear pattern that does emerge from 
our data is that the Co/Ni devices have stronger 1// 
noise than the NiFe devices, as illustrated by the rep- 
resentative spectra in Fig. [6J The white noise for the 
two device types is typically comparable (factor of 2.5 
difference in Fig. [6|, while the 1// component is typi- 
cally 10 times larger in the Co/Ni devices (factor of 15 
difference in Fig. IgI). Whether this systematic difference 
is due to the different materials in the STO free layer, 
the different magnetic anisotropics (which lead to differ- 
ent precessional trajectories), or to other factors is an 
important topic for future investigations. 

Several considerations rule out sources for the 1// 
noise other than the STO itself. Noise from the bias cur- 
rent source is filtered by the dc path of the bias tee, which 
has a bandwidth of 8 kHz, whereas we observed 1// noise 
up to much higher frequencies. Furthermore, an effect 
due to bias current noise would scale with Idvo/dli^l of 
the STO, but the devices shown in Fig. [6] follow the op- 
posite trend: the NiFe device has a larger duo/dI]j (by a 
factor of ~ 3) but weaker 1/ f noise than the Co/Ni de- 
vice. Noise from other sources such as the oscilloscope, 
preamplifier, or stray magnetic fields would affect all our 
measurements equally, which is not consistent with the 
reproducible differences among devices that we observed. 

Another trend emerges when we compare the measured 
SA linewidth, Ai^sA, with the value expected from the 
white noise level in Sj^{f)^ Ai^wh- (At this point we report 



Lorentzian linewidths for a measurement time of 1 s; see 
below for why it is important to specify both line shape 
and time scale.) Our NiFe devices have SA linewidths 
that are 1.1 to 1.3 times larger than expected from the 
white noise: Ai^sA = 11.2 MHz and Ai^wh = 10.4MHz 
for the device in Fig. p] (all values of Au here have an un- 
certainty of approximately ±0.5 MHz unless error bars 
on a plot indicate otherwise). In constrast, our Co/Ni 
devices have SA linewidths that are 2 to 3 times larger 
than expected from the white noise: Ai^sA = 10.0 MHz 
and Az^wh = 4.1MHz for the device in Fig. |6] As we 
describe next, this trend can be understood as a conse- 
quence of the different 1// noise strength in the different 
devices. 

Spectral line broadening due to 1// frequency noise 
is well known in the field of single-mode semiconduc- 
tor lasers. We first give a brief description of the key 
concepts and then apply them to our measurements in 
the next paragraph. White frequency noise in a laser, 
caused by spontaneous emission events, gives a spectral 
line having a Lorentzian shape [17] . Frequency noise hav- 
ing a 1// spectrum, caused by the fluctuating number 
of charge carriers in the semiconductor, gives an addi- 
tional contribution to the spectral line having a Gaussian 
shape [T8l[T9]. When the white and 1// contributions 
to the spectral line are comparable, the shape can be 
described by a convolution of Lorentzian and Gaussian 
profiles known as a Voigt function [I8l [19] . Importantly, 
the broadening depends not only on the white and 1// 
noise strengths, but also on the time scale of the mea- 
surement. Semiconductor laser spectra are typically mea- 
sured by interfering the light with a delayed copy of itself 
at a photodiode detector and modulating one arm of the 
interferometer (typically at 100 MHz) to avoid low fre- 
quency photodiode noise [20 . (The delay is achieved by 
placing a length of optical fiber in one arm of the interfer- 
ometer.) This method essentially translates spectral lines 
from optical frequencies to radio frequencies, where they 
can be measured with a conventional SA. The delay time 
Tdei in such a measurement sets a lower limit on the fre- 
quency noise: only Su{f) above / ^ 1/Tdei contributes to 
the spectral line [18]. Once T^ei is long enough that this 
lower limit lies in the 1// region of Sy{f)^ the spectral 
line will become broader as Tdei increases. A numerical 
study of cases where the white and 1// noise contribu- 
tions were comparable showed that the overall (Voigt) 
linewidth varies approximately logarithmically with de- 
lay time [T8] . 

To apply the semiconductor laser picture to an elec- 
tronic oscillator such as an STO, we must consider the 
appropriate time scale for an SA measurement. For a sin- 
gle SA sweep, the analog of Tdei is the sweep time over 
which the LO moves through a span around vq. The 
minimum sweep time for commercial instruments is typ- 
ically 1 ms. Comparing this time with the spectra in 
Fig. [6] we see that individual SA sweeps involve times 



over which the frequency noise in our STOs is not white. 
Moreover, several sweeps are usuahy averaged together 
and there is a dead time of ^ 100 ms between sweeps 
(required for restabiUzation of the LO). We typicahy av- 
eraged 10 sweeps to produce the final Sv{j^) that we fit 
to determine Az/sa, thus the total measurement time was 
^meas ^ 1 s. Certainly frequency noise below f ^ IRz 
cannot contribute to the spectral line, but the dead time 
prevents an exact mapping between T^ei in the laser mea- 
surement and Tmeas foi" our averaged SA measurements. 
We therefore proceed by developing a numerical model of 
a swept SA that simulates the signal processing occuring 
after the mixer, i.e., the steps that convert Vw{t) into a 
spectrum Sv{j^) averaged over multiple sweeps. 




Bin /_ Bin /+! *^ 



V(t) 




FIG. 7: (Color online) Illustration of how a swept SA gener- 
ates a spectrum from a voltage waveform. 

The process by which an SA converts a time domain 
signal into a frequency domain spectrum, illustrated in 
Fig. [7J includes various time scales that must be incor- 
porated into a model SA. For the ith frequency bin in 
the spectrum, the LO frequency of the mixer is fixed 
at z^LO during a segment time Tgeg (although SAs often 
sweep the LO continuously, here we consider it to be 
stepped discretely). The mixer IF output goes through 
a bandpass filter, whose width is set by the resolution 
bandwidth (RBW) of the SA, and then through a power 
detector [30]. The power detected during the ith seg- 
ment is the value of Sy for the ith frequency bin. For a 
spectrum spanning n frequency bins, the time to acquire 
a single spectrum is Tg^eep = ^^seg- When N spectra 
are averaged together, the time to restabilize z^lo at the 
first bin is typically much larger than Tgweep and thus 
the total measurement time T^eas is much larger than 
A/'Tgweep- Although this picture omits many details of 
the inner workings of an actual SA, it contains the rele- 
vant time scales that determine how non- white frequency 
noise affects the spectrum. 

We created a model SA that takes Fif(^) as the input, 
rather than V{t)^ since this is the data we recorded for 
our STOs. Thus rather than changing z^lo to generate 
different frequency bins, we changed the center frequency 



of the bandpass filter. After filtering each segment (using 
a filter with 3 dB points at half the RBW away from 
the center frequency and a rolloff of 36dB/decade), we 
computed the "power" in each bin by simply summing 
the squares of the values in the segment [31 . From each 
set of n segments we obtained a single spectrum, and 
we repeated the process to obtain A^ individual spectra 
for each value of Tgeg. We fit the individual Sv{^) to 
determine Az^sa for T^eas = ^^seg and we fit the average 
Sv{i^) to determine Az^sa for T^eas = A/'nTgeg (unlike 
in a real SA, there is no dead time between sweeps in 
our model SA). We used primarily the Voigt function 
to fit /Sy(z^) because it can fit lines that are Lorentzian, 
Gaussian, or any mixture of the two. We first present the 
Voigt results and then discuss fits using pure Lorentzian 
and Gaussian functions. 

Figure [S] shows Voigt linewidth vs. inverse measure- 
ment time for the Co/Ni device measured by (1) an ac- 
tual SA, and (2) our model SA applied to the same 50 ms 
Vw{t) waveform used to compute Sj,{f). For the model 
SA, we used n = 100 and varied T^eas by choosing values 
of Tgeg between 1 /is (shorter segments gave unreliable re- 
sults) and 500 /i5. We report both the mean from fits to 
10 individual spectra and the fit to the average of these 
10 spectra. For the actual SA, we varied T^eas by aver- 
aging with A^ = 1, 2, 5, 10, 100, and 1000, repeating the 
measurement five times for each value of T^eas hi order to 
report a mean value and estimate its uncertainty. The ac- 
tual and model SA results both show that Az^sa increases 
logarithmically with measurement time, and they agree 
quantitatively for T^eas = 1ms. When applied to the 
measured ViY{t) for the NiFe device, the model SA yields 
a weaker dependence of Az^sa on measurement time (the 
change in Az^sa barely exceeds the statistical uncertainty 
over the accessible range of Tmeas), which is consistent 
with the weaker 1// noise in Sj,{f) for this device. We 
also applied the model SA to numerically generated sig- 
nals. For signals having white frequency noise, Az^sA 
was independent of measurement time, whereas for sig- 
nals having Sy{f) comparable to that measured for our 
Co/Ni device, Az^sa '^s. 1/Tmeas had a slope similar to 
that seen in Fig. [8J Thus the model S A applied to signals 
having a range of Sy{f) from strictly white to strongly 
1// indicates that the dependence of linewidth on mea- 
surement time is a direct consequence of non- white fre- 
quency noise. 

As mentioned above, broadening due to 1// frequency 
noise is accompanied by a change in the shape of the 
spectral line [181 [IS ^ with the Voigt shape changing from 
mostly Lorentzian to mostly Gaussian as T^eas increases. 
The noise in our SA data for most values of Tmeas is too 
large to discern this trend clearly, i.e., a pure Lorentzian 
shape fits about as well as a Voigt shape. However, for 
the CoNi device at T^eas > 10 s a Voigt shape clearly 
fits the data better than a Lorentzian shape, as shown in 
Fig. [9J We also show a pure Gaussian fit to this line for 
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FIG. 8: (Color online) Dependence of Voigt linewidth on in- 
verse measurement time for the Co/Ni device. Model SA pa- 
rameters are RBW = 1 MHz, n=100,iV=10 (except for the 
longest segments where the 50 ms IF waveform was used up 
before reaching 10 spectra). Error bars for solid circles show 
the standard deviation of the mean for a set of N values with 
the same Tmeas- Actual SA parameters are RBW = 1 MHz, 
Tsweep = 1 ms, dead time between sweeps ^ 100 ms. Error 
bars for solid diamonds show the standard deviation of the 
mean for a set of 5 successive measurements with the same 
T 
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FIG. 9: (Color online) Spectral line for Co/Ni device with 
Tmeas = 10 s. Liucs show fits to Voigt, Gausslau, and 
Lorentzian functions. 



comparison, and we see that the Voigt function provides 
the best fit to the entire line. The Lorentzian shape gives 
a good fit to all spectra from the NiFe device, although 
we did not measure beyond T^eas = 1 s for the particular 
oscillation mode presented here. This is consistent with 
previous work on similar NiFe devices that did not find 
significant deviations from a Lorentzian line shape [2TJ 
[22]. 

In terms of the overall picture of STO spectral lines, 
our results mean there are two sources for a Gaussian 



contribution to line shape. The first involves the rate 
of relaxation to the stable precessional orbit, combined 
with the dependence of frequency on precessional am- 
plitude (the "frequency nonlinearity" intrinsic to STOs 
[7]), which sets a correlation time for phase fluctuations. 
When this correlation time is short compared to the ther- 
mal dephasing time the line shape is Lorentzian; in the 
opposite limit it is Gaussian [7 . In terms of frequency 
noise, this correlation suppresses Sy{f) at large / accord- 
ing to 



s.U) 



^wh 



1 + (2-/A)2 



(5) 



where rj is the relaxation rate ^ . The second source for a 
Gaussian line shape is 1// frequency noise, as described 
above, which can be understood as a correlation at long 
times that affects Sy{f) at small /. 

Distinguishing the two mechanisms for a non- 
Lorentzian STO line shape clearly requires more than a 
single SA measurement. Since the relaxation mechanism 
does not depend on measurement time, SA measurements 
over a wide range of T^eas could indicate whether the 1// 
mechanism is significant, but a quantitative conclusion 
about the relative contributions of the two mechanisms 
would require considerable care. Another approach is 
to use the autocorrelation function of the STO signal to 
measure the correlation at short times, as done in re- 
cent experimental work [S^ where the deviation from a 
Lorentzian line shape was less pronounced than that seen 
in Fig. [9J Finally, both the relaxation rate and 1// noise 
can be seen directly in Sjj{f) if it is measured over a large 
enough range of /. This approach has the advantage that 
each mechanism can be quantified separately. 

As mentioned above, in our measurements the lim- 
ited IF bandwidth suppressed Sy{f) above / ^ 20 MHz. 
While this prevents a direct measurement of 7^, it does 
set a lower bound of 7^/27r '^ 20 MHz (see Eq. ([5|). From 
the theory of the relaxation mechanism ^ i9^ , the condi- 
tion for a negligible Gaussian contribution can be written 
as rj/ 27: ^ Ai^wh- Since the measured values of Ai^wh 
(10.4 MHz for the NiFe device; 4.1MHz for the Co/Ni 
device) are less than 20 MHz, we can conclude that the 
line shape for both devices would be close to Lorentzian 
in the absence of 1// noise. The strongly non-Lorentzian 
line in Fig. ^ for the Co/Ni device can be unambiguously 
attributed to 1// frequency noise. 

Although non-white frequency noise has not been di- 
rectly measured in previous STO experiments, recent re- 
ports indicate it is probably an important effect in de- 
vices beyond the two types considered here. The line 
"jitter" and SA linewidth increasing with measurement 
time reported in MgO nanopillar STOs [24] can both be 
interpreted as evidence for 1// frequency noise, and the 
authors suggest a possible mechanism for such noise that 
is specific to their particular devices. In other nanopillar 
STOs containing MgO [25 or metallic [26 barriers, the 



DFT of selected short segments yielded linewidths much 
smaller than those found from either DFT or SA mea- 
surements averaged over long times. Beyond these pub- 
lished reports, we and others have noticed when watch- 
ing the non-averaged SA display that some devices show 
larger trace-to-trace jumps in center frequency than oth- 
ers, an observation that may be explained by varying 
amounts of 1// noise in the devices. 



CONCLUSION 

We measured the output of two types of STOs in the 
time domain, using the IF ouput of an SA to access sig- 
nals well above 10 GHz. We presented two techniques for 
obtaining the power spectrum of frequency noise, Su{f)^ 
and showed the advantage of combining them to yield an 
averaged Sjy{f) over a wide range of /. The 1// noise 
we observed indicates that theoretical models based on 
thermal noise sources are insufficient for our devices over 
times longer than about 1 /is. We also measured spec- 
tral linewidths in our devices using both an actual SA 
and a numerical model that allowed shorter measurement 
times. For the devices with the strongest 1// noise, we 
found that SA linewidth increased, and the line shape be- 
came significantly non-Lorentzian, as measurement time 
increased. These results imply that SA measurements of 
STOs should be accompanied by measurement time val- 
ues so that (1) comparisons can be made among various 
STOs measured by different researchers, and (2) non- 
Lorentzian line shapes can be correctly interpreted. Al- 
though the consequences of our results for STO appli- 
cations depend on the time scales involved, we expect 
that measurements of Sjy{f) will allow more accurate 
predictions of performance than SA measurements alone 
for many applications. The detailed picture of oscillator 
noise provided by Su{f) will also help in distinguishing 
among different physical origins for the noise. 
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